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Abstract. Hot Jupiters are new laboratories for the physics of giant planet 
atmospheres. Subject to unusual forcing conditions, the circulation regime on 
these planets may be unlike anything known in the Solar System. Characterizing 
the atmospheric circulation of hot Jupiters is necessary for reliable interpretation 
of the multifaceted data currently being collected on these planets. We discuss 
several fundamental concepts of atmospheric dynamics that are likely central to 
obtaining a solid understanding of these fascinating atmospheres. A particular 
effort is made to compare the various modeling approaches employed so far to 
address this challenging problem. 



1. Introduction 

An exploding body of observations promises to unveil the meteorology of giant 
planets around other stars. Because of their high temperatures, short orbital 
periods, and likelihood of transiting their stars, the hot Jupiters are yielding 
their secrets most easily, and we now have constraints on radii, composition, 
albedo, dayside temperature structure , and even day-night temperature dis- 



Harrington et al.l 



tributions for a variety of planets (e.g., Knutson et al. 20071 : Cowan et al. 20071: 
' l l2006l . l2007l : ICharbonneau et al.ll2002l . l2005l . l2007l : lDeming et al 



2005,, ,20061 ). A knowledge of atmospheric dynamics is crucial for interpreting 



these observations. 

Understanding the atmospheric circulation of any planet is a difficult task, 
and hot Jupiters are no exception. The difficulty results from the nonlinearity of 
fluid motion and from the complex interaction between radiation, fluid flow, and 
cloud microphysics. Turbulence, convection, atmospheric waves, vortices, and 
jet streams can all interact in a complex manner across a range of temporal and 
spatial scales. Furthermore, radiative transfer and dynamics are coupled and 
cannot be understood in isolation. For example, the equator-to-pole tempera- 
ture contrasts on terrestrial planets depend not only on the rate of latitudinal 
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energy transport by the circulation but also on the way the atmospheric radi- 
ation field is affected by the advected temperature field. Cloud microphysics 
complicates the problem even more, because the large-scale radiative proper- 
ties of clouds depend on the microphysics of the cloud particles (particle num- 
ber density, shape, size distribution, and absorption/scattering properties). We 
thus have a non-linear, coupled radiation-hydrodynamics problem potentially 
involving interactions over 14 orders of magnitude, from the cloud-particle scale 
(~ 1 ^m) to the global scale (~ 10^ m for hot Jupiters). 

It is worth emphasizing here the important difference between the com- 
plexity of planetary atmospheres and the relative simplicity of stars, exempli- 
fied by the main sequence in the HR diagram. While specifying a few global 
parameters — such as mass, composition, and age — is typically sufficient to un- 
derstand the key observable properties of stars, this is generally not the case 
for planetary atmospheres. It is possible that a significant observable diversity 
exists even among a group of extrasolar planets which share similar global at- 
tributes. Understanding such complexity is a challenge for extrasolar planetary 
science. 

A proven method for dealing with the complexity of planetary atmospheres 
is the concept of a model hierarchy. Even a hypothetical computer model that 
included all relevant processes and perfectly simulated an atmosphere would not, 
by itself, guarantee an understanding of the atmospheric behavior any more than 
would the observations of the real atmosphere (jSeld 2005). This is because, in 
complex numerical simulations, it is often unclear how and why the simulation 
produces a specific behavior. To discern which processes cause which outcomes, 
it is important to compare models with a range of complexities (in which various 
processes are turned on or off) to build a hierarchical understanding. 

For example, east- west jet streams occur in the atmospheres of all the plan- 
ets in our Solar System. The study of these jets involves a wide range of models, 
all of which have important lessons to teach. The most idealized are pure 2D 
models that investigate jet formation in the simplest possible conte xt of hori 



zontal, 2D turbulence int e racting with planetary rotation fe.g.. iWilliami 



Yoden and Yamadal Il993l : ICho and Polvanil Il996bl : E uang and RobinsonI 



1978; 



1995 



Sukorianskv et al. 20071 ). Despite the fact that these simplified models exclude 



vertical structure, thermodynamics, radiation, and clouds and provide only a 
crude parameterization of turbulent stirring, they produce jets with several sim- 
ilarities to those on the planets. Next are one-layer "shallow-water-type" models 
that allow the fluid thickness to vary, which introduces buoyancy waves, alters 
the vortex interaction le ng ths, and hence changes the details of je t formation 
(ICho and Polvanil ll996al lbr iScott and Polvanil l2007l : IShowmanI 120071 ). 3D mod- 



els with simplified forcing allow the investigation of jet vertical structure, the 
interaction of heat tra nsport and jet f o rmation, and the 3p stability of jets 
to various instabilities (ICho et al.l l200ll : IWilliamsl Il979l . l2003l : ISchneided l2006l : 
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Lian and Showmanll2007l ). Such models suggest, for example, that Jupiter and 



Saturn's superrotatin£| equatorial jets may require 3D dynamics. Finally are 
full 3D general-circulation models (GCMs) that include realistic radiative trans- 



^Superrotation is defined simply as a positive (eastward) longitudinally averaged wind velocity 
at the equator, so that the atmospheric gas rotates faster than the planet. 
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fer, representations of clouds, and other effects necessary for detailed predictions 
and comparisons with observations of specific planets. The comparison of simple 
models with more complex models provides insights not easily obtainable from 
one type of model alone. This lesson applies equally to hot Jupiters: a hierarchy 
of models ranging from simple to complex will be necessary to build a robust 
understanding. 

Here we review our current understanding of atmospheric circulation on hot 
Jupiters. We first describe basic aspects of relevant theory from atmospheric 
dynamics; this is followed by a detailed comparison of the equations, forcing 
methods, and results obtained by the different groups attempting to model the 
atmospheric circulation on these fascinating objects. 

2. Basic considerations 

Despite the complexity of atmospheric circulation, there exists a 40-year history 
of work in atmospheric dynamics on solar-system planets that can guide our 
understanding of hot Jupiters. 

Differences between hot Jupiters and solar-system giants: Jupiter, Saturn, Uranus, 
and Neptune rotate rapidly, with rotation periods ranging from 10 hours (Jupiter) 
to 17 hours (Uranus). This implies rotationally dominated flows. In contrast, 
hot Jupiters are ex pected to rotate synchronously with their orbital periods 
( Guillot et al.|[l996l ). at least when orbital eccentricity is small, implying rota- 
tion periods of 1-5 days for the known transiting planets. Some hot Jupiters 
will therefore be rotationally dominated, while on others the Coriolis forces will 
have more modest importance. Note, however, that even "slowly" rotating hot 
Jupiters rotate much faster than Titan and Venus (rotation periods of 16 and 
243 Earth days, respectively). On most hot Jupiters, rotation will have a major 
impact in modifying the flow geometry. 

The forcing on hot Jupiters differs greatly from that on solar-system gi- 
ants. First is simply the amplitude of the radiative forcing: HD209458b absorbs 
and reradiates almost 260, 000 W m^^ on a global average, as compared to 
240 Wm~^ for Earth, 14 Wm^^ for Jupiter, and 0.7Wm~^ for Uranus and 
Neptune. A second difference is that, on the giant planets in our Solar Sys- 
tem, the flow seems to be driven by convection, thunderstorms, or baroclinic 
instabilities with small length scales (see below). Jupiter, Saturn, and Neptune 
have comparable instrinsic and total heat fluxes, allowing interior convection to 
greatly affect the atmosphere; the convective zone extends to the visible cloud 
layers at ~ 1 bar on these planets. For hot Jupiters, however, a statically stable 
radiative zone extends to pressures of 100-1000 bars, separating the observable 
atmosphere from the convective interior; furthermore, the interior convective 
forcing is lO'^ times weaker than the stellar insolation. Thus, the hemispheric- 
scale stellar forcing is likely to be more dominant on hot Jupiters than for the 
giant planets in our Solar System. Finally, the expected radiative times above 
the photospheres on hot Jupiters are extremely short — days or less — whereas 
those on solar-system giants are years to decades. All these differences will 
impact the flow in ways that remain to be understood. 

Despite our ability to list these differences in forcing, the way the circulation 
responds to these differences is often subtle. The planets in our Solar System 
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demonstrate that the relationship between forcing and response is nontrivial. For 
example, Neptune and Uranus exhibit similar wind patterns but their forcings 
differ greatly. Neptune's intrinsic heat flux is 1.6 times the absorbed sunlight, 
implying a strong role for interior convection; on Uranus, in contrast, the intrin- 
sic flux is < 10% of the total flux, implying that external forcing is predominant. 
Moreover, the small (30°) obliquity on Neptune implies the existence of a regu- 
lar day-night pattern where the equator receives more sunlight than the poles; 
in contrast, the 98° obliquity on Uranus means that the poles receive greater 
sunlight than the equator, and that much of the planet experiences decades of 
near-continuous sunlight followed by decades of darkness. This Uranus-Neptune 
comparison suggests that a rotating-stratified atmospheric flow can respond to 
forcing in a manner which is not simply dictated by the forcing geometry. 



Flow length scales: Theoretical work has demonstrated the importance of two 
fundamental length scales in atmospheric dynamics. First, the Rhines length, 
defined as (JJ/P)^^'^, is the scale at which planetary rotation causes east-west 
elongation (jets). Here U is a characteristic horizontal wind speed and /3 = 
2^}cos(p/a is the latitudinal gradient of planetary rotation (Q is the planetary 
angular rotation rate, (f) is the latitude, and a is the planetary radius). Tur- 
bulence at small scales tends to be horizontally isotropic, but at length scales 
approaching the Rhines scale, turbulent structures grow preferentially in the 
east-west direction rather than the north-south direction. This typically leads 
to a flow exhibiting a banded structure with east-west jet streams. 

Second, the Rosshy deformation radius, defined as NH/f, is the scale at 
which pressure perturbations are resisted by the Coriolis forces (/ = 20, sin (j) 
is the Coriolis parameter, H is the scale height, and N is the Brunt- Vaisala 
frequency, i.e., the oscillation frequency for a vertically displaced air parcel in 
a statically stable atmosphere). Vortices often have horizontal sizes near the 
deformation radius. And baroclinic instabilities — a form of sloping convection 
that can occur in the presence of horizontal temperature contrasts, converting 
potential energy to kinetic energy — typically generate turbulent eddies with a 
horizontal size comparable to the deformation radius. Such instabilities are the 
primary cause of midlatitude weather on Earth and Mars and may be i mportant 
in dr i ving t he jet streams on Jupiter, Saturn, Uranus, and Neptu ne ( William^ 
19791 . 120031 : iFriedson and Ingersol]|[l987l : iLian and ShowmanI [20071 ^. 

On Jupiter and Saturn, the deformation radius and Rhines length are, re- 
spectively, ~ 2000 and ~ 10,000km (i.e., much less than the planetary radius). 
This explains why these planets have narrow jets, numerous cloud bands, and 
small vortices. On most planets, including Earth's atmosphere, Earth's oceans, 
and the atmospheres on Mars, Jupiter, and Saturn, there exists a large popula- 
tion of flow structures with sizes near the deformation radius. 

In contrast, most hot Jupiters have modest rotation rates, large scale heights 
(up to ten times that on solar-system giants as a result of the high temperature), 
and strong vertical stratification (a result of the dominant external irradiation) . 
This means that, on most hot Jupiters, the Rhines length and Rossby deforma- 
tion radius are close to the planetary radius. Unlike Jupiter and Saturn, the 
dominant flow structures ori hot Jupiters should therefore be planetary in scale 
(jShowman and Guillotl[200l : iMenou et al"ll2003l ^. 
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The length scales discussed here will strongly influence the atmospheric 
energetics. In a 3D turbulent fluid, vortex stretching occurs readily, generating 
small-scale flow structures from large-scale ones and causing energy to cascade 
downscale. In a quasi-2D fluid such as an atmosphere, however, vortex stretching 
is inhibited, and energy instead experienc es an inverse c a scade that t ransfers 
energy from small to large length scales (jPedloskvl [TqSTI : IVallid 120061 ') . (The 
simplest example is vortex merger, which gradually produces larger and larger 
vortices over time.) The extent to which the inverse cascade occurs depends 
on the length scale of energy injection relative to the planetary size: if energy 
injection occurs at small scales (as may be the case on Jupiter and Saturn), the 
inverse cascade occurs readily. On the other hand, if energy injection occurs at 
planetary scales (as may occur on hot Jupiters), the inverse cascade could be 
less dominant (but see §3.2 below). 



Wind speeds: Despite progress, there is still no general theory for what sets 



the wind speeds on the giant planets in our Solar System. The speeds on all 
four giants exceed those on Earth, which presumably relates to the lack of a 
surface (and its associated friction) on the giant planets. More puzzling is the 
observation that the wind speeds on Neptune exceed those on Jupiter despite 
the fact that Neptune's total atmospheric heat flux — which ultimately drives 
the circulation — is only 4% of that on Jupiter. A lesson is that the mean speed 
involves an equilibrium between forcing and damping (e.g., friction). Fast winds 
can occur with weak forcing if the damping is also weak; alternatively, slow 
winds can occur in the presence of strong forcing if the damping is strong. The 
damping mechanisms on solar-system giants (which include turbulent mixing 
and radiative cooling to space) seem to be relatively weak, leading to fast winds. 
However, it is not fully understood how these mechanisms translate into the 
hot Jupiter context, particularly since the strength of the damping (as well as 
the forcing) depends nonlinearly on the flow itself. Careful theoretical work can 
shed light on how damping and forcing interact to produce an equilibrated flow 
speed relevant for hot Jupiters. 

Despite the above caveats, the known dynamical link between winds and 
temperatures can be invoked to make useful statements about the winds. In 
rotating atmospheres, horizontal temperature contrasts a re linked to vertical 
differences in horizontal wind via the thermal-wind relation ( Holton 20041 ) . (This 
is often defined in the context of geostrophic balance but it can be extended to 
a more general balance.) On terrestrial planets, the speed is zero at the surface. 
So, given an estimate of the equator-to-pole temperature difference, the thermal- 
wind relation can provide a crude estimate of the actual wind speeds. The forcing 
and damping processes are relatively well understood for terrestrial planets, and 
fully nonlinear 3D circulation models do a reasonably good job of simulating the 
structure and speeds of the winds. 

For giant planets, one can profitably apply the thermal- wind relation to 
obtain the difference in horizontal wind speed between (say) the cloud deck and 
the base of the weather layer (thought to be at ~ 10 bars for Jupiter) (e.g., 
IngersoU and Cuzzil Il969l ). The difficulty on giant planets is that, unlike the 



terrestrial planets, the speed at the base of this layer is not necessarily zero 
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the planet may exhibit a so-called "barotropicQ' wind that penetrates through 
the convective molecular envelope. In this case, the total wind at the cloud 
tops would be the sum of the thermal-wind (or "baroclinic" component) and 
the barotropic wind. 

In the context of hot Jupiters, several numerical simulations have been 

performed that are driven by the intense dayside heating and nightside cooling 

( Showman and Guillotll2002 : Cooper and Showmanll2005l l2006l : lDobbs-Dixon and Lin 



20071 : iLangton and LaughlinI 



20071 ^. These simulations are effectively constrain- 



ing the thermal- wind component of the flow (the height-variable component that 
exists in the radiative zone above ~ 100-1000 bars). All of these simulations 
obtain wind speeds of ~ Ikmsec"^, in agreement with order-of-ma gnitude esti- 



mates using the thermal-wind relation (jShowman and Guillotll2002i ). All of the 



published simulations assume that the deep barotropic wind component is zero. 
It is worth bearing in mind, however, that the barotropic wind component is 
unknown, and if it is strong then the actual wind speed could differ in magnitude 
and geometry from that suggested by these simulations. Detailed comparisons 
of simulations with observed lightcurves and spectra of hot Jupiters can provide 
constraints on the extent to which the models are capturing the correct flow 
regime. 

Temperatures: Insight into temperature patterns may be gained with timescale 
arguments. Suppose Trad is the radiative timescale (i.e., the time for radiation 
to induce large changes in the entropy) and Tadvect is the timescale for air to 
advect across a hemisphere (e.g., time for air to travel from dayside to nightside 
or equator to pole). If r^ad ^ ''"advect! then large temperature contrasts are 
expected, whereas if Trad ^ '^advect, then temperatures should be horizontally 
homogenized. 

To illustrate the usefulness of these arguments. Table 1 lists these timescales 
and temperature contrasts for the tropospheres of the terrestrial planets. The 
numbers show general agreement with the above arguments. The radiative time 
constants are long on Venus, intermediate on Earth, and short on Mars. Re- 
garding transport, we must distinguish longitude from latitude. On Earth and 
Mars, the planetary rotation dominates the transport of air parcels in longitude, 
so parcels cycle between day and night with a timescale of 1 day. Venus rotates 
slowly, but winds and rotation nevertheless transport air from dayside to night- 
side over a period of days to months. The numbers in the Table imply that 
"^rad S> Tadvect,lon for Vcnus and Earth; in agreement, longitudinal temperature 
variations are small on these two planets. (They reach ~ lOK on Earth, but 
this results largely from the uneven continent-ocean distribution; a continent-free 
"aquaplanet" Earth would exhibit even smaller longitudinal temperature varia- 
tions.) Latitudinal transport timescales (Table 1) imply that Trad > ''"advect.iat for 
Venus, Trad ~ Tadvectjat for Earth, and r^ad < Tadvect.lat for Mars. In qualitative 



barotropic fluid is one where the contours of constant pressure and density align, while 
a baroclinic fluid is one where these contours misalign. The convective interior has nearly 
constant entropy and is hence barotropic. Under sufficiently rapid rotation this leads to the 
Taylor-Proudman theorem in the convective interior, which states that east-west wind is con- 
stant in columns that penetrate the molecular envelope in the direction parallel to the rotation 
axis. This is called a barotropic wind. 
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agreement with expectations, equator-to-pole temperature contrasts are small 
on Venus, intermediate on Earth, and large on Mars (Table 1). 



Table 1. Timescales and temperature contrasts on terrestrial planets 



Planet 


Trad 


Tadvcct,lon 


ATlon 


'^advect.lat 


ATiat 


Venus 


years 


days 


<fewK 


months? 


fewK 


Earth 


weeks-months 


1 day 


- lOK 


weeks 


20-40 K 


Mars 


days 


1 day 


~ 20-40 K 


days-weeks 


~ 50-80 K 



Values refer to tropospheres (0.005 bars for Mars, 1 bar for Earth, and ~ 10 bars on 
Venus). In this table, "day" refers to Earth day (86400 sec). Refs: iGierasch et al.l 
(flQQTl ): IPeixoto and Oor^ (fl99l ; iRead and Lewii (pOOl . 



What about giant planets? The four solar-system giant s exhibit latitu- 
dinal temperature contrasts above the clouds of only ~ 5K (jingerso 111990). 
This results from several processes. First, in the interiors, convective mixing 
times are short (~decades) but t he radiative timescales ar e geological, so the 
interior entropy is homogenized ( Ingersoll and Pored 19781 ). On Jupiter, Sat- 
urn, and Neptune (which have internal convective heat fluxes comparable to 
the absorbed solar flux), the top of the constant-entropy interior is close to the 
photosphere; thus, we expect nearly constant photospheric temperatures across 
the planet. Second, on all four giants, the radiative timescales above the clouds 
are long (years to decades), and so even relatively sluggish circulations are ca- 
pable of homogenizing the temperature patterns. The latter mechanism may be 
particularly important on Uranus, which lacks a significant interior heat flux. 

Hot Jupiters differ from solar-system giants because the radiative timescales 



den below a radiative zone extending to ~ 100 -1000 bars 
Guillot an d Showman"'2002': 'Bodenhei mer et al.ll200ll . l2003l : 



Chabrier et al, 2004; Iro et al. 2005 ). 



(Guillot et al. 


1996; 


Baraffe et al. 


2003; 


re expected photo- 



sphere pressure of ~ 0.01-1 bar. Large day-night temperature differences are 
therefore possible at the photosphere depending on Tj-ad and Tadvect there. Cur- 
rent estimates, which derive from radiative calculations that neglect circulation, 
suggest that r^ad is extremely pre ssure dependen t, ranging from ~ 1 hour at 

Considering wind speeds 
"^advect ^t pressures of 



20051 ). 



10"'^ bars to ~ 1 year at 10 bars (jlro et al. 
ranging from 0.3-3 km sec~^ then implies that r^d 
~ 0.1-1 bar. Above this level, day-night temperature differences should be 
large, while below this leve l the temperature differences could become small 
(|Showman and GuillotI [20021 1. 

A major caveat, of course, is that one must have knowledge of the circu- 
lation to estimate Tadvect- Even r^ad can depend on the circulation through its 
dependence on temperature and composition (which in turn depend on the cir- 
culation). In the absence of other information, the above arguments are thus 
generally insufficient to make rigorous predictions of temperature contrasts on 
planets, since independent timescales are assigned to radiation and hydrodynam- 
ics while the two processes are in fact intimately coupled, as mentioned earlier. 
Nevertheless, the arguments can be used a posteriori to understand the temper- 
ature patterns occurring in observations — or in detailed nonlinear numerical 
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simulations of the circulation — when these observations and/or simulations are 
sufficient to define Tadvect and Trad- 



Models of atmospheric circulation on hot Jupiters: 
and differences 



Similarities 



A variety of approaches have been used to study the atmospheric circulation 
on hot Jupiters. These include solutioi is of the global two-diinensio nal shallow- 
water equations on a rotating sphere (jLangton and Laughlin 2007), the global 
three-dimensional primitiv e equation^ on a rotating sphere ( Showman and Guillot 
2002l : ICooper and Showma n 2005, 2006), the global equiv alent-barotropic formu- 
lation of these primitive equations on the rotating sphere ( Cho et al. 20031 . 20071 ). 
the two-dimensio nal Navier-Stokes eq uations neglecting rotation in an equato- 
rial cross section ( Burkert et al. 20051 ) . and the three-dim ensional Navier-Stokes 
equa tions omitting the polar regions of a rotating sphere (jPobbs-Dixon and Lin 
20071 ). The forcing in all of these models is simplified; no models yet include 
realistic radiative transfer. 



3.1. Work of Showman, Cooper and collaborators 

IShowman and GuillotI ( 20021 ) . [Cooper and ShowmanI ( 20051 ) . and lCooper and Showman 
(j2006 ) solved the 3D primitive equations on the whole sphere for parameters rel- 
evant to HD209458b. The 3D approach is motivated by the expectatio n that the 



radiative time constant varies by orders of magnitude in the vertical (jiro et al 



|2005| ) , which in the presence of the day-night heating gradient would cause pat- 
terns of temperature and wind that vary both vertically and horizontally in an 
inherently 3D manner. 

The primitive equations are the relevant equations for large-scale flow in 
statically stable atmospheres whose horizontal dimensions greatly exceed the 
vertical dimensions. This is expected to be true on hot Jupiters, which have 
horizontal scales of 10''-10'^ m but atmospheric scale heights of only 200-500 km 
(depending on temperature and gravity), leading to a horizontahvertical aspect 
ratio of 20-500. This large aspect ratio allows the vertical momentum equation 
to be replaced with local hydrostatic balance, meaning that the local vertical 
pressure gradient dp/dz is balanced by the local fluid weight pg, where p, p, g, 
and z are pre ssure, density, gravity, and height (see for example the derivation in 



Holtonll2004l . and §3.5 below). Lateral and height variations in these quantities 



can still occur and evolve dynamically with the flow. The primitive equations 
admit the full range of balanced motions, gravity (i.e. buoyancy) waves, rota- 
tionally modified (e.g., Rossby and Kelvin) waves, and horizontally propagating 

sou nd waves, but they filter out v ertically propagat ing souiid waves. ^ 

, [Sh owman and GuillotI (j2002l ) . [Cooper and ShowmanI (j2005l ). and lCooper and ShowmanI 

(|2006l ) : breed the flow using a Newtonian heating/cooling scheme, which param- 
eterizes the thermodynamic heating rate as (Tgq — r)/ri.ad, where Teq is the 
specified radiative-equilibrium temperature profile (hot on the dayside, cold on 



■^The leading-order dynamical equations so lved in meteorology and climate studies are called 
the primitive equations (e.g.. lHoltonll20d^ : see §3.1 and §3.5. 
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the nightside), T is the actual temperature, and Trad is the radiative-equihbrium 
timescale ( a function of pr essure). The vertical structure of T^q and Trad were 
taken from [iro et al. ( 20051 ): the day- night difference Tgq was a free parameter 
that was varied from 100 to 1000 K. This formulation, of course, does not repre- 
sent a rigorous treatment of radiative transfer, but it provides a physically mo- 
tivated means to investigate how the expected day-night thermal forcing would 
drive the atmospheric circulation. 
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Figure 1. Simulated temp erature (greyscale) and horizo ntal winds (arrows) 
for three pressure levels from lCooper and ShowmanI (|2006l ). Planetary param- 
eters of HD209458b were adopted. 



Figure 1 shows the temperature (greyscale) and winds (arrows) for three 
layers (10, 100, and 1000 mbar from top to bottom, respectively) after a simu- 
lated time of 1000 Earth days for a case with an assumed day-night difference in 
Teq of 1000 K. The imposed heating contrast leads to winds of several kmsec"^, 
and the day-night temperature contrasts reach 500 K or more. The flow struc- 
ture has a strongly 3D character. At the top (pressures less than ~ 100 mbar), 
the radiative time constant is shorter than the time for winds to advect the 
flow, and the temperature pattern tracks the heating — hot on the dayside and 
cold on the nightside — despite the fast winds. At greater pressure, however, 
the radiative time constant is longer, and dynamics blows the hottest region 
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Wavelength (fim) 




Figure 2. Theoretical infrared spectra (l eft) and lightcu r ves in Spitzer IRAC 
and MIPS ba nds (right) as calculated bvlFortnev et alJ ()2006( ) from the 3D 
simulations of lCooper and ShowmanI ()2006 ) for HD209458b. The dayside and 
nightside spectra (labelled "secondary eclipse" and "transit" , respectively) dif- 
fer greatly due to the distinct vertical temperature-pressure structures in each 
region. Case 0, 1, and 2 refer to different assumptions about the chemistry. 
In the lightcurves, note the ph ase shift of the peak flux from the time near 
secondary eclipse (at 180°). See lFortnev et al.l (|2006f ). 



of the atmosphere downwind from the substellar point by 20-60° of longitude 
depending on altitude. In agreement with the Rhines-length and deformation- 
radius arguments presented in §2, the simulations exhibit only a small num- 
ber of broad jets. Interestingly, the development of a su perrotating (eastward) 
equat or ial jet is a robust featu r e in a ll the simulations of IShowman and Guillot 
( 2OO2I ) ■ ICooper and ShowmanI ( 20051 ). and ICooper and Showman ( 20061 ). The 
mechanism for driving the jet seems to involve equatorward pumping of eddy 
momentum that occurs in conjunction with the longitudinally varying heating. 
This process has p reviously been observed in simplified two -layer models from 
the Earth sciences (jSuarez and Duff v 199^ ISaravananI Il993l ) . 

By pushing the temperature pattern far from radiative equilibrium, the cir- 
culation has major implications for the infrared lig htcurves and spectra . Using 

(I2OO6I) calcu- 



Fortnev et al 



a two-stream plane-parallel radiative-transfer code, 
l ated i nfrare d lightcurves from the 3D temperature patterns of lCooper and Showman 
( 2OO5I . [20061 ) . These calculations predict that the peak IR emission leads the sec- 
ondary eclipse by ~ 2-3 hours in most Spitzer IRAC bands (Fig. 2). This results 
from advection of the temperature field by the dynamics, which displaces the 
hottest regions to the east of the substellar point (Fig. 1). Intriguingly, the 
predicted offset in the flux peak is similar to that in the observed lightcurve for 
HD189733b, but the simulations overpredict the day-night flux difference and 



fail t o produce to cold region to the west of the antistellar point (jKnutson et al 
20071 ). 



The dynamics also drives the vertical temperature profile far from radiative 
equilibrium. As air advects from nightside to dayside in these simulations, the air 
aloft warms faster than air at depth (because of shorter radiative time constants 
aloft), producing a quasi- isothermal region on the dayside despite the absence 
of such a feature in pure radiative-equilibrium models. This leads to dayside 
and nightside spectra (Fig. 2, left) that differ substantially from the predictions 
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of ID radiative-equilibrium models (iBurrows et al 
Seager et al1l2005l : iFortnev et alJ[200i 



20051 : iBarman et aD l2005l : 



3.2. Work of Cho, Menou and collaborators 

In geophysical fluid dynamics, an important distinction is made between flows 
that are fundamentally 3D and flows whose essence is largely 2D (and hence 
can be captured in a 2D model), even though they are of course occurring in a 
3D atmosphere. For example, the Earth's midlatitude troposphere experiences 
baroclinic instability, which is a fundamentally 3D process and requires a 3D 
model to capture. Large parts of the Earth's stratosphere (the region overlying 
the troposphere), on the other hand, lack baroclinic instabilities because of the 
strong static stability. Many aspects of the stratosphere — such as the dynamics 
of the polar vortex, Rossby wave breaking, and turbule nt mixing — have been 
mode led to great a.dvantage with 2D or quasi-2D models ( Juckes and Mclntvrd 
Il987l : [jucked 19891 : iPolvani et al.lll995l ). Because of their strong static stabilities 
and larg e Rossby d e forma t ion ra dii (whi ch tend to discourage baroclinic insta- 
bilities), ICho et al.1 (|2003l . l2007l ^ and M enou et al.l (j2C)03) have proposed that 
the basic flow regime of hot Jupiter atmospheres may likewise be successfully 
captured in a quasi-2D, one-layer formulation. 






--HJO -201} 2U0 400 

Eastward Velocity |ni/s| 



Figure 3. [Left] Equatorial and polar views of potential vorticity (a flow 
tracer) in a specific hot Jupiter model from Cho et al. (2003, 2007). Note the 
prominent circumpolar vortices formed as a result of potential vorticity con- 
servation. [Right] Corresponding zonally averaged wind profile, characterized 
by a small number of broad jets (three in this case) . 



Adopting this point of view, ICho et al.l (|2003l . 120071 ') solved the equivalent- 
barotropic formulation of the primitive equations on the rotating sphere. These 
equations are a one-layer version of the primitive equations and can formally 
be obtained by integrating the primitive equations vertically. The vertical in- 
tegration leads to equations for the column-integrated horizontal velocity, layer 
thickness, and mean temperature as a function of longitude, latitude, and time. 
The advantage of this approach is that, by reducing the vertical resolution. 
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high horizontal resolution can be achieved, allowing numerical solutions to re- 
solve small-scale turbulent interactions that are difficult to capture in a full 3D 
model. These equations (and the formally related shallow-water equations) are 
well-suited to the study of two-dimensional jets, vortices, and horizontally prop- 
agating gravity and Rossby waves. These simplified equations have a venerable 
history of successful use in studying d ynamical processes in the plan etary and 
atmospheric sciences (for a review see IVasavada and Showman |2005| ) . On the 
other hand, the full three-dimensional nature of the flow is not represented. 
These equations exclude 3D processes such as baroclinic instability, vortex tilt- 
ing, and vertical wave propagation. 



R ather than exciting the flow only with a day-night heating contrast. ICho et al 



( 20031 . 20071 ) also initialize their simulations with small-scale turbulent stirring. 
Their flows subsequently evolve to states containing meandering polar vortices, 
breaking Rossby waves, and turbulent mixing (Fig. 3). One may question the 
relevance of small-scale initialization for hot Jupiter atmospheres, since they 
are intensely heated on large, hemispheric scales. In considering this question, 
however, one should remember that coupling between small and large scales in 
nonlinear flow is a delicate, and still not completely understood, issue. As an 
example, consider once again Uranus: it is dominantly forced by insolation on a 
hemispheric scale and yet the rotationally constrained flow is clearly not a sim- 
ple reflection of this external forcing. It is significant that equivalent-barotropic 
simulations with small-scale initial turbule nce are able to r eproduce the domi- 
nant features of Uranus' atmospheric flow ( Cho et 311120071 ). Such simulations 



demonstrate that regardless of the type and scale of forcing, rotation plays a 
dominant role in shaping the large-scale flows on planets. 



2.01 r 




Mm 1300 1400 1EDD K 2 4 6 8 

orbital phase 



Figure 4. [Left] Equatorial views of four successive temperature maps for 
the same model as shown in Fig. 3. Large-scale moving atmospheric fea- 
tures, in particular circumpolar vortices, may lead to significant observable 
variabihty. [R ight] Exan i ples of thermal infrared phase curves according to 
the models of ICho et al.l (|2007D . From top to bottom, curves for adopted 
mean wind speeds of 800, 400, 200 and 100 m/s are shown. At higher wind 
speeds, contributions from coherent moving atmospheric features become ap- 
parent, shifting and d eforming the phase curves away from a regular shape 
(|Rauscher et all 120081 ). 
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It is imp ortant to realize t hat t he flow is stirred only at the start of the 
simulations in Cho et al. ( 20031 . [2^0071 ) . The resulting flow in these dynamically 
initialized studies can be extremely time-variable, especially when the mean 
winds in the turbulent initial condition are strong. If such a behavior occurs 
on a hot Jupiter, then substantial variability in the secondary ecli pse depths 



of a transiting hot Jupiter may be expected ( Rauscher et al.ir2007bl ). In addi- 



tion, signatures of large-scale moving atmospheric features may b ecome apparent 



(Fig. 4) in the thermal infrared phase curves of these objects (jRauscher et al 
20081'). In agreement with Rhines-length and deformation-rad ius arguments 



(jShowman and Guillotll2002l: ICho et alJl2003l : iMeriou et aDl2003l1 and numerical 
simulations by ISh owman and GuillotI (|2002l ) and ICoooer and Showma^ (l2005l . 
20061 ). the flow usually consists of a small number of broad jets (Fig. 3). 



3.3. Other work 

iLangton and LaughlinI ( 20071 ) performed global numerical simulations of the at- 
mospheric flow on hot Jupiters using the traditional shallow- water equations, 
which govern the dynamics of a thin layer of hydrostatically balanced, constant- 
density fluid on a rotating sphere. Mathematically, these equations can be 
viewed as a special case of the equivalent-barotropic equations in isentropic 
coordinates when R/cp = 1, where .R is the gas constant and Cp is the specific 



heat at constant pressure ( Cho et al] l2007l ). However, because the density in 
the shallow-water equations is constant (whereas it is dependent on pressure in 
the equivalent barotropic equations), th e temperature i s undefined and there is 
no therr nodynamic energy equation ( see |Pedloskvlll987l . Chapter 3). Note that 
Eq. 1 of iLangton and LaughlinI (|2007l ). which they write as an energy equation, 
is actually the shallow- water mass continuity equation for the layer thickness: 



dh 
'dt 



+ V • (/iv) 



(1) 



where h is layer thickness (a representation of pressure), t is time, and v i s 
horizontal velocity of the shallow- water layer. ILangton and Laughlin ( 20071 ) 
force Eq. 1 by including a Newtonian-relaxation source term (/leq — ^)/7"rad) where 
/leq is the equilibrium thickness (thick on dayside, thin on nightside), which 
adds mass on the dayside and removes it on the nightside. This is a reasonable 
zeroth-order me thod to model the effects of radiat ive cooling on the shallow- 
water flow (e.g.. [Polvani et ah 19951 : Showman 20071 ). A rigorous representation 
of thermodynamic forcing and temperature requires the use of a more general 
set of equations (such as the equivalent-barotropic, primitive, or compressible 
Navier-Stokes equations) that, unlike the shallow-water equations, have distinct 
energy equations. In the planetary context, the shallow-water layer could be 
interpreted as a thin, stable "weather l ayer" , possibly overlying a deeper and 
denser interior. Langton and Laughlin ( 20071 ) adopt a radiative time constant 
of ^8 hours, intended to re present conditions expected near the photosphere 
according to Iro et al. ( 20051 ). These forced simulations are initialized from rest|l| 



ILangton and LaughlinI (j2007l ) also pr esent unfor c ed sim ulations initialized with small-scale tur- 
b ulence, simi l ar to the simulations of lCho et al.l l|2003f) . In this case, they obtain a flow similar 
to ICho et al.) (1200^ 1. 
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Langton and LaughlinI ( 20071 ) find that, when the obhquity is assumed zero, 



their forced flows quickly reach a steady state exhibiting a westward equato- 
rial jet with speeds approaching Ikmsec"^ and day-night differences in layer 
thickness reaching the imposed equilibrium value of ~ 40%. The height field 
becomes distorted from the equilibrium value, /leq, so that the thickest and 
thinnest regions are offset from the substellar and antistellar points. Interest- 
ingly, despite the development of a westward equatorial jet, these thin and thick 
regions are displaced to the east rather than to the west of the antistellar and 
substellar points, respectively ; the reasons for this phenomenon remain unclear. 
Langton and LaughlinI ( 20071 ) also performed simulations with a high obliquity, 



wherein the starlight alternately heats one hemisphere and then the other; the 
dynamics here becomes periodic with the orbital period, as expected. In all 
these simulations, the flow length scales are broad, similar to the previous stud- 
ies discussed earlier. 

Several authors have also performed investigations in limited-area domains. 
Burkert et al. ( 20051 ) perform a 2D study (with coordinates longitude and height) 



using the Navier-Stokes equations without rotation in a section of an equatorial 
plane. They parameterize the heating with a flux-limited diffusion scheme using 
Rosseland mean opacities and an imposed temperature distribution at the top 
of the model typically ranging from 100 K on the nightside to ~1200K at the 
substellar point. This scheme induces heating on the dayside and cooling on the 
nightside. Their standard domain extends 180° in longitude and ~6000 km verti- 
cally with free-slip wall boundary conditions. Because of the side wall boundary 
conditions, the flow cannot develop planet-encircling jet streams (as occurs for 
exarnple in the simula, t ions of[Sh owman and G uillot 2002; Cooper and ShowmanI 
20051 : ICho etalll2003l . l2007l : [Pobbs-Dixon and Linll2007i l. Instead, an overturn- 



ing circulation develops with rising motion near one wall, sinking motion near 
the other, and lateral currents that connect the two. The nominal opacity used 
for calculating the diffusion coefficients in the radiative scheme is based on that 
expected for grains in the interstellar medium, although simulations were also 
performed with opacities ranging from 10^^ to lO'^ times the nominal values. 
The strength of the flow depends on the assumed opacity; in the nominal cases, 
the wind speeds reach several kmsec"^ and day-night temperature differences 
reac h ^700 K at the expec t ed ph otosphere pressure. 



Dobbs-Dixon and LinI ( 20071 ) investigated the atmospheric circulation by 



solving the 3D Navier-Stokes equations in the low- a nd mid-latitud e s of a rotat- 
ing sphere (the polar regions were omitted). As in iBurkert et al. ( 20051 ). they 



parameterized the radiative-transfer with flux-limited diffusion using Rosseland- 
mean opacities. In addition to the Coriolis force, they included the centrifugal 
force in the equation of motion; this differs from the approach taken in most 
global atmospheric circulation models, which typically account for the gravita- 
tional relaxation of the planetary interior to t he planetary rotation by absorb- 
in g the centrifugal force i nto th e gr avity (see Holt od l2004l. p. 12- 13). Similar 



to lShowman and GuillotI (|2002l ') and [CooEerandShowman (2005. l2006l ). their 
simulations develop an eastward jet at the equator and westward jets in midlat- 
itudes. The vertical thickne ss of these jets varies substantially with longitude. 



Dobbs-Dixon and LinI (120071) flnd t hat the flow strength depends on the assumed 



opacitv. as in iBurkert et al. ( 20051 ); in their nominal cases, wind speeds reach 
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several km sec ^ and day-night temperature differences at the expected photo- 
sphere reach ~ 700-800 K. 

The studies described s o far all focus on flow near the photosphere. In con- 
trast, [Koskine^^er^l] ( 20071 ) performed 3D simulations of the thermosphere and 
ionosphere of hot Jupiters using the 3D primitive equations. In their simulations, 
the domain extended from ~ 2 /ibar to ~ 0.04 nbar or less. The simulations in- 
cluded molecular heat conduction and diffusion, which are important due to the 
large molecular mean- free paths at the sub-/ibar pressures investigated in this 
study. Simplified schemes for photochemistry, stellar extreme ultraviolet (EUV) 
heating, and radiative cooling by Hj^ were also implemented. The emphasis was 
on planets from 0.2-1 AU from their stars (in contrast to the other investiga- 
tions described in this review, which tend to emphasize planets at < 0.1 AU). 
The simulations developed strong flows of ~1 km see" ^ from the substellar point 
to the antistellar point; temperatures often reached 10, 000 K or more depend- 
ing on assumptions about the cooling. The wind, temperature, and chemistry 
patterns have implications for the rate of atmospheric escape to space. 



3.4. Speculations on the causes of model differences 



Published simulations of atmospheric circulation on hot Jupiters have simi- 
larities as well as differences. Despite the range of approaches, all the pub- 
lished studies agree that the flow on typical hot Jupiters will exhibit only a 
small number o f broad jets, consistent with Rhines scale and deformation ra- 
dius argum ents ("Showman and Guillot"2002'; 'Cooper and Showman"2005', '200^. 



Menou et al. 2003; Cho ct al. 20il3,_2007; Langton and Laughlin 2007; Dobbs-D ixon and Lin 
2007i ;i. As discussed earlier, this agreement results from the modest rotation 



rates, high temperatures, and strong static stabilities of hot- Jupiter atmospheres. 
(Indeed, this bodes well for characterizing the dynamical state with spectra and 
lightcurves, which sample globally or near- globally.) Published simulations that 
attempt to predict the wind speeds also generally agree that the baroclinic com- 
ponent of the flow (i.e., the difference in the flow speeds between the photo- 
sphere and the base of the radiative zone) may reach ~ 1 kin sec~^ or more 
jShowman and Guillot' '2002'; 'Coope r and Show man' '2005', '200^; iBurkert et aD 
[2005 : .Langton and Laughhn, ,2007 : Dobbs-Dixon and Lin 2007 ). 

On the other hand, the models also exhib i t a number of important differ- 
ences. Characteristic features of the Cho et al. ( 2003 . 20071 ) simulations include 
development of large polar vortices, extensive turbulent mixing, and time vari- 
ability, if the global average wind speed is large and/or the da y-night hemispheric 
tempe rature difference is weak. In con trast, the simulations of Showman and GuillotI 
(120021 1. ICooper and ShowmanI (l2005l . l2006l ). and iLaneton and LaughlinI (|2007l ) 
exhibit steadier flow beha vior and te nd to lack the high degree of turbulent 
mixing that occurs in the Cho et ahl simulations. While differences in model 
resolution and numerical methods surely contribute, one cause of these differ- 
en ces seems to be the is sue of whether the initial flow contains turbulence (as 



m 



Cho et al 



... ^x..,2 003. 2007*), i n addition to a global-scale day- n i ght forcing (as in 

Showma n and Gui llot 200z; 'Coop er and Showman|[2005l . 20061 : Langton and Laughlin 
2007; Do bbs-Dixon and Lin 2007). Certainly, if the initial condition contains 
strong, closely-packed vortices (as in ICho et all ) , the subsequent process of vor- 



tex mergers and planetary wave radiation naturally induces turbulent mixing 
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and produces, as an end s tate, a small number of large vortices that typically 
end up near the poles (see Cho and Polvanilll996b l). On the other hand, when 
the energy source comprises global-scale day-night forcing, the dominant length 
scales are large (~planetary radius), leading to a global-scale flow. 

Therefore, a key distinction between the various models published so far has 
to do with the problem of understanding how the stellar forcing is expressed in 
the atmospheric flow itself. As we have already commented, this is a subtle issue, 
even for solar-system giant planets. However, there are good reasons to believe 
that progress on this issue is possible in the future. First, on the modeling side, 
detailed schemes to treat the radiative transfer in hot Jupiter atmospheres are 
being developed, so that more reliable treatments of the non-linear radiation- 
hydrodynamics involved in expressing this forcing can be implemented. These 
modeling efforts may thus be able to shed some light as to what the nature of the 
flow forcing in hot Jupiter atmospheres is. Secondly, observations may provide 
important guidance for the modelers, as is the case even in solar-system atmo- 
spheric studies|£| This is particularly true of the issue of atmospheric turbulence. 
All solar-system atmospheres are turbulent. For hot Jupiters, there is an expec- 
tation that the combination of turbulence and large fundamental atmospheric 
scales may lead to significant variability in various observables. Observations 
may thus help distinguish a rather steady atmospheric behavior from a more 
variable one, and in doing so establish a distinction between the two modeling 
approaches (and associated forcing schemes) discussed above. As has always 
been the case with solar-system studies, progress in understanding hot Jupiters 
is likely to involve repeated iterations between observations and modeling. 

Another difference in published stu dies is the development of robust equato- 

rial superrotation in the 3D simulations ( Showman and Guillotll200"^ : ICooper and Showman 
2005 . 20061 : Dobbs-Dixon and Lin|[2007l ). whereas the one-layer calcul ations pro- 



duce both westward (i.e., subro tating) and eastward equatorial flow (jCho et al, 



20031. l2007l: iLangton and Laug hlin 200^). Regarding the superro t ation , the s tud- 
les bylSh owman and Guillotf(|2002i ) and ICooper and ShowmanI (|2005l . |2006| ) in- 
volved a range of different assumptions for the initial temperature profile, the 
day-night equilibrium temperature difference, and the radiative time constants. 
Some simulations were initialized with a strong westward jet to determine if 
it could be retained in the simulation. Intriguingly, all of these simulations 
produced a strong eastward jet at the equator. This suggests that the super- 
rotation is a robust phenomenon, at least within the context of the adopted 
input parameters and forcing approach. More study is needed to understand 
this phenomenon. 

The forcing in all published models i s simplified; no mode ls ye t include real- 

i stic r epresentations of radiative transfer. Burkert et al. ( 20051 ) and lDobbs-Dixon and Lin 
( 20071 ) adopted a diffusion approach, which is accurate in the deep, optically 
thick atmosphere (and hence is commonly used in the interior portion of evolu- 



^Large-scale jet streams can exhibit shear instabilities that generate small-scale turbulence, but 
this phenomenon has not been obviously apparent in the radiatively forced studies performed 
to date. 

''Even advanced GCMs, which are treated as much as possible as first-principle models, generally 
include numerous free parameters which are adjusted to match observed atmospheres. 
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tion models); however, the radiative transfer is non-diffusive above the photo- 
sphere, so this approach loses accuracy in the observable atmosp here (i.e. at and 

above the photosphere) . In contrast, the shallow- water st udies oflLangton and Laughlin 
2007 ) and t he 3D studies of lShowman and GuillotI ( 20021 ) and lCooper and Showman 



2005 . I2OO6I ) adopted a Newtonian heating/cooling scheme. This scheme has 



enjoyed a long history of successful use in process studies of planetary climate. 
However, it effectively involves a linearization of the heating rate around the 
radiative-equilibrium state and neglects the fact that the radiative equilibrium 
temperature and timescale can depend on the atmosphere's dynamical response, 
particularly when actual ter nperatures are far from the prescribed radiative- 
equilibrium values. Finally, ICho et all ( 20031 . |2007| ) did not include radiative 
forcing at all but instead performed adiabatic simulations driven by a combi- 
nation of initial turbulence and large-scale pressure deflections. This approach 
allows an exploration of how the flow behavior depends on the energy (which is a 
tunable input parameter) and therefore provides a useful baseline. On the other 
hand, there is no ability to predict mean flow speeds, and the absence of radia- 
tive forcing means that these simulations cannot capture aspects of the dynamics 
that crucially depend on such forcing. The adopted approaches, while simpli- 
fied, nevertheless all serve the intended purpose of providing plausible sources 
of energy injection. Inclusion of realistic radiative transfer is an important goal 
for future modeling work. 

Interestingly, infrared lightcurve observations could allow an observational 
constraint on the wind patterns, at least for planets with a photosphere deep 
enough for the radiative time to be comparable to the advection time. To vary- 
ing degrees, the simulations performed to date suggest that dynamics can dis- 
tort the temperature pattern in a complex manner, so that the temperature 
field does not simply track the stellar heating pattern ("Showman and Guillqt 




,Cooper and ShowmanI 12005 . 2006.: Cho et al. 2003 . 2007: Burkert et al 
Dobbs-Dixon and Lin l2007l : iLangton and Laughhnl l2007l ) . For example. 



the hottest and coldest regions may be displaced from the substellar and antis- 
tellar points, respectively, leading to phase shifts in infrared lightcurves relative 
to that expected without dynamics. While there will likely be degeneracies 
in interpretation (i.e., multiple flow patterns may explain a given lightcurve), 
such data will nevertheless provide powerful contraints to distinguish among the 
various models. 



3.5. The validity of hydrostatic balance 

A scaling analysis demonstrates that local hydrostatic balance is approximately 
valid for the large-scale flow on hot Jupiters. It is important to emphasize 
that the local-hydrostatic-balance assumption in the primitive equations derives 
from the assumption of large aspect ratio and not from any assumption on wind 
speed. However, the fact that estimated wind speeds on hot Jupiters are several 
km see" ^, which is close to the 3-kmsec~^ speed of sound in these atmospheres, 
suggests that we consider the validity of the primitive equations in the hot- 
Jupiter context. The full Navier-Stokes vertical momentum equation can be 
written 

dw „ 1 dp ^ , 

— + vVw = — TT- - 5 + cos (2) 
ot p oz 



18 



Showman, Menou, and Cho 



where w is vertical wind speed, t is time, v is the horizontal wind velocity, and u 
is the east-west wind speed. The background static hydrostatic balance is irrel- 
evant to atmospheric circulation and can be removed from the equation. Define 
p = po{z) + p' and p = po{z) + p' where po and po are the time-independent 
basic-state pressure and density and, by construction, dpo/dz = —pog. Primed 
quantities are the deviations from this basic state caused by dynamics. Substi- 
tuting these expressions into Eq. [2l we can rewrite the equation as 

dw „ 1 / dp' , \ ^ , 

— + V • = — - pg] + 2un cos (/) (3) 



dt p \dz 

where the basic-state hydrostatic balanced has been subtracted off. The terms in 
parentheses give only the flow-induced contributions to vertical pressure gradient 
and weight (they go to zero in a static atmosphere). 

For local hydrostatic balance to be a reasonable approximation, the terms 
in parenthesis must be much larger than the other terms in the equation. The 
magnitude of dw/dt is approximately w/t, where r is a flow evolution timescale, 
and the magnitude of v • Vw is the greater of uw/L and nP' / H where L is the 
horizontal flow scale and H is the vertical flow scale. For global-scale flows, L ~ 
10^-10*^ m and r ~ 10 ^-10^ sec. The large- scale flo w varies vertically over a scale 
height H 300 km ( Showman and Guillot 200_^; ICooper and ShowmanI |2005| . 



20061 : iDobbs-Dixon and Linll2007l 1. For the simulated flow regime in Cooper and 



Showman (2005, 2006) and Dobbs-Dixon and Lin (2007), v ~ u ~ 3kmsec~^, 
g ~ lOmsec"^, $7 ~ 2 x lO^^sec"^, and w ~ 10-100 msec^^. With these 
values, we find that dw/dt < 10~^msec~^, uw/L < 0.03msec~^, w'^/H < 
0.03msec~^, and Vtu ~ O.lmsec"^. In comparison, for a hot Jupiter with day- 
night temperature differences of several hundred K, the flow-induced hydrostatic 
terms p'g/ p and p^^dp' /dz are each ^ lOmsec^^. 

This analysis implies that, for global-scale hot-Jupiter flows, the greatest de- 
parture from hydrostaticity results from the vertical Coriolis force, which causes 
a ~ 1% deviation from hydrostatic balance. The acceleration terms on the left 
side of Eq. [3] (which are necessary for vertically propagating sound waves) cause 
a ~ 0.3% deviation from hydrostatic balance. Hydrostatic balance is thus a rea- 
sonable approximation for the large-scale flow. Nonhydrostatic effects of course 
become important at small scales, and it is conceivable that these effects interact 
with the large-scale flow in nontrivial ways. For hot Jupiters, the acceleration 
terms on the left side of Eq.[3]only become important for structures with vertical 
or horizontal scales less than ~ 30 km and ~ 500-1000 km, respectively. In a nu- 
merical model that solved the full Navier-Stokes equations, the grid resolution 
would have to be substantially finer than these values for the nonhydrostatic 
behavior to be accurately represented. A full Navier-Stokes solution with a 
coarse resolution would effectively be resolving just the global-scale hydrostatic 
component of the flow. 

3.6. The need for model validation 

Due to the complexity of processes involved in modeling the coupled radiation- 
hydrodynamics of a planetary atmosphere, it is customary in the fleld of geo- 
physical fluid dynamics and atmospheric sciences to implement robust validation 
schemes for any simulation tool. Typically, both radiative transfer modules and 
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dynamical cores (i.e., "primitive equation solvers") are subjected to standard 
test cases or inter-comparisons aimed at quantifying their performance and ac- 
curacy in conserving important quantities, such as energy and enstrophy ("vor- 
tical energy"). For example, a classical benchmark test for dynamical cores is 
the Held-Suarez test, which isolates baroclinic dynamics from radiative transfer 
issues by using a simple, Newtonian cooling scenario (Held and Suarez 1994). 

The issue of model validation is a particularly critical one for hot Jupiter 
studies. Earth and planetary scientists have access to detailed information on 
the atmospheres they study, allowing serious model flaws to be rapidly iden- 
tified. This is not the case with hot Jupiter atmospheres. While information 
available from observation of these distant atmospheres is rapidly increasing, it 
still remains subject to degeneracie s in interpretation. This is so even if "eclips e 
maps" eventually become available (IWilhams et al.ll2006l : lRauscher et al.ll2007al ). 
As a result, rigorous model validations using test cases and Solar System exam- 
ples will be critical. Successful model validations can be achieved by explic- 
itly reproducing known features in the Solar System planet atmospheres (e.g. 
ICho et al.ll2007li. using simulation tools which ha.ve been well-tested previously 
fe.g. IShowman and GuillotI l2002l : ICho et al.l[2003 : ICooper and ShowmanI l2005l ) 
or by designing appropriate inter-comparison tests tailored for hot Jupiter prob- 
lems. 



4. Summary 

The study of hot Jupiter atmospheres is maturing. While data of increasing 
quality are being collected, atmospheric models are also being refined to help 
build a robust, hierarchical understanding of these unusual atmospheres. Indeed, 
one of the main motivations for studying these atmospheres, which is also a 
main source of difficulty, is the unusual physical regime that characterizes them. 
Hence, hot Jupiters represent new laboratories for studying the complex physics 
of giant planet atmospheres. In this way, they offer the promise of extending 
the boundary of comparative planetology well beyond the solar-system planets. 
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